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Orbital Stability for the Schrodinger Operator 
Involving Inverse Square Potential 

Georgios P. Trachanas, Athens* 

Nikolaos B. Zographopoulos, Athens ^ 


Abstract 

In this paper we prove the existence of orbitally stable standing waves for the critical 
Schrodinger operator, involving potential of the form The approach, being 

purely variational, is based on the precompactness of any minimizing sequence with respect 
to the associated energy. Moreover, we discuss the case of the presence of a Hardy energy 
term, in conjunction with the behavior of the standing waves at the singularity. 


1 Introduction 

We are concerned with the concept of orbital stability of the standing wave solutions for the 
following semilinear Schrodinger equation 

(1.1) ipt + Ap + (^ ^2 1^ ^ x€R^, 

where p : x M —>■ C, N > 3 and 2 < q < + 2. 

A standing wave is a solution of the form 

'4>{x,t) = A G M, u : R^ —>■ M, 

therefore, is reduced to the stationary equation 

—^ yA + An — |nk“^n = 0, x G M^. 

/ \x\ 

It is well-known that Hardy’s inequality 

[ \^dx < [ iVnkdx, for all n G iV > 3, 

RN kl Jm 

is closely related to elliptic and parabolic equations involving inverse square potentials. The 
optimal constant c* := { ^ 2 ^ ) is the natural borderline separating existence from nonexistence. 
We know from |35j that the linear heat equation with an inverse square potential, 

Q 

(1.3) Ut = Au + - — prU, 

\x\ 
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admits a global solution for c < c*, both in the bounded and the whole space case. On the other 
hand, there are no solutions for c > c*, even locally in time, due to instantaneous blow-up. In 
the critical case c = c*, a global solution is defined in [331 EH 135] > but the functional framework 
is more complicated. 

There is a quite extent literature related to the Schrddinger operator equipped with an inverse 
square potential. He := —A — c|x|“^. Regarding the stationary problem, we shall refer to some 
papers related to the semilinear equation 

(1.4) HcU = f{x,u). 

In |18] . the author exhibits some existence, uniqueness and regularity results on various types 
of solutions for the case 0 < c < c*, posed on a bounded domain containing the origin, where 
f{x, u) = u'^ + tg{x), 5 > 0 is a smooth, bounded function and t > 0. The case t = 0 on a ball of 
appears in [5]. 

In the strictly subcritical case c G (0, c*), we refer indicatively to |3T] on some existence and 
nonexistence results in concerning the case f{x,u) = rt^-2. We also mention [28j . 

JV+2 

where the case f{x,u) = K{x)u^-^ is considered, yielding a nonexistence result for c > c*, and 
K{x) is a positive and bounded weight. In [T7|, some existence results are obtained in 

4 

where f{x,u) = a{x)u -|- \u\’^-^u + g{x,u) for certain functions o, g. In the context of kernel 
estimates for the semigroup , 0 < c < c*, we refer to see also |2] for the subcritical case 
of the potential. Analogue bounds on the Schrddinger heat kernel for the case of the Laplace- 
Beltrami operator on some Riemmanian manifold can be found in m- A result on the long-time 
behavior for the magnetic Schrddinger operator with the critical potential is obtained in |8]. 
Finally, we refer to [T] for the spectral analysis, in of the operator 

La vU := —Art — — riu, 

\x\^ 

where 0 < < 1 and rj satisfies certain growth conditions. We also mention m for some 

estimates on the moments of the negative eigenvalues of the Schrddinger operator in the critical 
case. 

In |29j , the existence of radially symmetric ground state solutions is established for the following 
semilinear equation 

(1.5) —/S.U + V{\x\)u = Q{\x\)u^~^, a; G 

where F(|a;|) is nonnegative and presents a singular behavior at the origin of order G 

[2, A). Analogue results are obtained in [30], by the same authors, for an equation similar to 
m, driven by the p-Laplacian. In both cases, the procedure is based on proving compact 
inclusions for certain weighted Sobolev spaces. Also, we refer to the variational approach of m 
for the case R(|x|) > —(c* — a)\x\ a > 0, R(|x|)|xp -|-oo, whenever |x| —>■ 0 or \x\ —oo. 

Finally, we refer to |3] for the case of a nonnegative potential V(x) = where a > 0 

and I G Z. 

In the time-dependent case, nonlinear equations of the general form 

(1.6) iV't-k A'i/'-k R(x)V’+/(x,-0) = 0, 

arise in various contexts of mathematical physics (see for example [26] and the references therein). 
Perhaps, in many circumstances, the orbital stability of the standing wave solutions is the only 
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extracted information on the asymptotic behavior of nonlinear phenomena governed by (jl.Op . 
When V{x) = 0, we refer to the stability result in [11] (see also [27]). In the same paper, a 
Hartree type equation is considered, equipped with a potential of the form V(x) = , 

where ci,...,Cm are positive constants, and xi,...,Xm are given poles. The case V G 
can be found in H [26]. We especially point out [22] where, among other, the existence and 
stability is obtained for a singular potential V(x) ~ |x| 6 G ( 0 , 2 ), near zero. In [THl 123] ( 2 ^. 

the authors prove the orbital stability for weighted nonlinearities of certain singularity, that is, 
they consider the case /{x^'ij)) = Q{x)\'il}\‘^~^'il}, where Q is allowed to behave like G (0,2). 

Concerning singularities at infinity, we refer to the instability results in HU [I3] for the case 
of harmonic potentials V{x) = |xp. A case of quasilinear equations is treated in |24j . In the 
subcritical case I4(|x|) = c|x|“^, c < c*, we refer to [38], where the long-time behavior is studied 
by applying scattering theory in the fI^(M^)-setting. On some Strichartz type estimates for the 
linear Schrodinger and the linear wave equation, respectively, equipped with a subcritical inverse 
square potential, we refer to | 6 ]. 

As it turns out, as far as we know, the results obtained here are new and they cover the critical 
inverse square potential, involving the best constant, on the whole space. 

In Section 121 we formulate the functional framework. The proper norm for the setting of problems 
m and () 1 . 2 I) is the following 

(1.7) ||u||^ := hm(/Bc(u) - A£(u)) + ||u||^ 2 (KAr), 

where 

Ib^{u):= [ \Vu\^dx-(^^) [ 

J\ ^ J Jb^ pI 

is the Hardy functional on the complement of the ball of radius e and centered at the origin, and 
the surface integral 

Ae{u) := ^ f I'UpdS', 

^ J\x\=e 

represents the Hardy energy at the singularity. This unexpected energy formulation was hrst 
introduced in the paper [33] in order to overcome a functional difficulty emerged by the eigen- 
solutions presenting the maximal singularity. Further applications to linear singular parabolic 
problems of the above consideration can be found in [321154] . 

The present approach is based on the existence of critical points for the following functional 

E{u) := lim(lBc(u) - A^{u)) [ |u|^dx, 

2 elo ^ q Jrjv 

on certain subsets of the energy space H with prescribed L^-norm. Namely, we define the En¬ 
sphere of H, 



for a given 7 > 0 , and we study the following minimization problem 


( 1 . 8 ) 


u G r, 

E{u) = min {E{z) : z G T, z is radial} , 
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that is, we are interested for radial ground state solutions. The set of all global minimizers of 
(II.8p is denoted by 

:= {u € T : u is a solution of (|1.8p } . 


In Section [3l we recall the notion of orbital stability and we prove the main result which is stated 
in the following: 

Theorem 1.1 Assume that 2 < q < -^ + 2, N > 3, and j > 0 is a given constant. Then, the 
set S.y is orbitally stable. 

Based on a proper transformation as in [33p [34] , we derive the behavior of the standing waves at 
the origin. More precisely, we establish that the minimizers of (11.81) behave exactly as 
at the singularity, thus they do not belong to In this case, we address the appearance 

of the Hardy singularity term at the origin. 

Theorem 1.2 For each 7 > 0, every minimizer of \1.^) behaves at the origin like . 


At this point we mention that the first nonexistence in result was given in |21j . concern¬ 

ing the corresponding to (11.21) linear problem. The first exact description of the behavior at 
the singularity, and thus nonexistence in the classical Sobolev spaces, was given in [331 ED for 
all minimizers. In almost all cases, this behavior corresponds to the appearance of the Hardy 
singularity term as a well-defined real number (positive or negative depending on the case); the 
Hardy functional is also well-dehned and finite as a principal value (also positive or negative de¬ 
pending on the case). Exception to this is the case of the so called /c-improved Hardy functional 
(see m), where both the associated Hardy functional and the Hardy singularity term tend to 
infinity. Finally, we mention that results concerning the behavior at the singularity, by a different 
method, may be found in [20] . 

Furthermore, we confirm the orbital stability for the following Schrodinger equation 


(1.9) 


•11-4 A I ^ ~ ^ 1 r \ 

i\x\ tctArc-b ( —^— I ^^ = -h[x)w, 


X G 


for some h decaying to zero at infinity, which also presents a Hardy-type energy. The somehow 
unexpected fact is that this energy appears in an additive way to the total energy, as an effect 
that comes from infinity, and maybe represents the main part of it. The result is based on the 
arguments of [33l[34]- If we consider the corresponding to (II.9p linear heat equation, 


( 1 . 10 ) 


1 1—4 
\x\ Ut 


X G 

M^, t > 0, 

u{x, 0) 

= ^^ 0 ( 3 ;), 

for 

X G 

u{x, t) 

—^ 0, 

for 

|x —>■ 00 , t > 0 


then, using similarity variables, we may prove the existence of the Hardy-type energy. This 
energy comes from infinity, is additive to the total energy and constitutes the main part of it. 

In Sectional we extend Theorem 1 1.1 1 bv removing the hypothesis of radial symmetry on xf and u. 
Nevertheless, in the general case, the cost we have to pay is the introduction of a certain weight 
function on the nonlinearity. Note that there is no possibility for a ” non-weighted” Hardy-Sobolev 
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inequality to hold, thus it is unclear to us if standard methods (e.g. see [9]) may be applied. On 
the other hand, in Subsection KT\ we prove a weighted (or improved) Hardy-Sobolev inequality. 

More precisely, we prove the orbital stability of the following equation 

~= 0, 

on X M, for weight functions g : —>■ [0, +oo) decaying, with a certain rate, to zero at 

inhnity. 


2 Space Setting 


In this section, we introduce a more convenient function space in order to slide over the effect of 
the Hardy term c*|x|“^. We define the space Ti to be the completion of the C'“(M'^)-functions 
under the norm 


( 2 . 1 ) 


|2 

m 


V = 


/ \x\ dx + 


|x| 


Notice that, setting 

_ N—2 

(2.2) u{x) = l^l ^ 'y(x). 


in ()2.1I) . the first integral simplifies the Hardy functional, since 





\SIu\^dx 


N-2y 


rdx. 


at least for C'“(R^)-functions. In the sequel, we denote transformation (12.21) by u = T(u), which 
is an isometry between the spaces X = L^(M'^) and X = L^(M'^,d/i), dg = \x\~^^~'^^dx. Our 
setting is based on this equivalence, as we shall see below. 

It is clear that 7i, equipped with the scalar product 

(2.3) (ui, U 2 )'H := Re / VuiVu 2 (ix + Re / \x\~^^~‘^^viV 2 dx, 

Js.^ Jrn 

is a well-defined real Hilbert space. An equivalent definition of H is the following: 


Lemma 2.1 Every function v, such that HuH-^ < oo, belongs to %. 


Proof. We adapt the arguments of [33] in our case. First, we prove the result for v € 
For £ > 0 small enough and t > 0, we define the cutoff function Pe(t) S C'o(M+\{0}), 0 
as: 


Pe:(t) 


' 0, 0 < t < £^, 

(log 1/e)-^ log(t/£^), £^ <t <£, 
< 1, e < t < i, 

(loge)-Mog(te2), 

. 0 , t>^. 


<Pe<l, 


For a fixed v G ||u ll"^ < oo, we define v^(x) = ps{\x\)v{x). Note that 


iVpedxDI = ^, Ce := (log(l/e)) \ 

\x\ 
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for X G Ae U , where = {x G < \x\ < e} and = {x G : e ^ < |x| < e ^}, 

e e 

being zero otherwise. Then, we have 

WVe-vWu < 2 / \V Pe{\x\)f Ivl"^ dx 

JAeUAi 

e 

(2.4) +2 [ \x\-^^-'^'> {I-Pef {\Vv\^+ \v\^)dx. 

£ 

We will prove that letting e J, 0, the integrals in (12.4|) tend to zero. The only one that is delicate 
is the hrst. We have 

f |x|-W-2) |Vps(|x|)p |xpdx < C||x||^ / c^—= C||x||^(log(l/e))"\ 

J As J ^ 

and ^ 

I iVpedxDI^ |xpdx < C'llxll^ / c^—= C||x||^(log(l/e))“\ 

■JAi j e~^ ^ 

£ 

and these tend to zero as e —>■ 0. 


Finally, we will prove that for any v with ||x||^ < 00 , v is approximated by a sequence of 
bounded functions in %. Indeed, if we define Vn as follows 


„ _ / ^(^)> if l^(^)l ^ 

"" \ n, if |u(x)| > n, 

then, we have 

(2.5) Ibnll'H = [ |Vxp dx + f |x|“i^“^i |up dx < 00 . 

Jc„ JCr, 

Now, it is clear that the sets Cn = {x G : |u(x)| > n} form a monotone family and the 
measure tends to zero in the limit n ^ 00. This means that the above integral goes to zero as 
n ^ 00 and the proof is complete. ■ 

Now, the approximation of the functions v^, which vanish around the origin, by C'“(M'^\{0}) is 
standard. 


Corollary 2.1 The set C'“(M^\{0}) is dense in TL. 

We introduce the space H to be the isometric space of TL under T, that is, H is dehned as the 
completion of the set 

{u = |xr^x : V G C'o°°(M^)} = r(Co°°(M^)), 

under the norm 

(2.6) N‘^{u) := [ |x|“i^“^i|V(|x|^“rt)pdx + f |npdx. 

Jr^ Jrn 

Then, the following is immediate: 

Corollary 2.2 The set C^{R^\{0}) is dense in H. 

Note also that 

(2.7) ||x||^ = ||u||^ , for all u G C^(M^\{0}). 

Arguing as in [33], the norm of H is related to the Hardy functional by means of the formula 

dEZD. 
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3 The Radial Case 

3.1 Space Properties 

In the following, 2* := 2N/{N — 2) stands for the critical Sobolev exponent. 

Lemma 3.1 Let Tir be the subspace ofLL consisting of radial functions. Then, 

TLr ^ |x| ~^dx), 

with compact inclusion for any 2 < q < 2*, N > 3. 

Proof. Let Vn{r), r = |a:|, be a bounded sequence of C'“(M'^\{0})- functions in PL] without loss 
of generality we assume that Vn converges weakly to 0. We will prove that 

f N—2 f N—2 f N—2 

(3.1) / \x\ ^'^\Vn\'^dx= / \x\ ’^^\Vn\^dx+ / |x| \Vn\^ dx ^ D, 

Jm.^ Jbi Jbi 

as n —)> oo. We claim that, for Vn in PL{Bi), it holds Vn S H^{Bi C M^). This is clear from the 
following 

lbn||«(Bi) ~ / r\v'J^dr+l r\vn\^ dr. 

Jo Jo 

Using now the compact imbedding 


(3.2) 

we obtain 


H^{Bi C M^) ^ L^{Bi C M^), for any p £ (1, +oo). 


/ \x\ \vn\^ dx ^ 0, for any p€(l,+oo), 

Jb, 


as n ^ oo. Then, for p large enough. 


L 


x\ 2 \Vn\'^ dx = / 

Bi JBi 


\X\ 2p ||3;| p y I? ylx 


\x 


^ ( / kl 2 (p-g) 

\Jb-^ j \jb^ 

The first integral in the right hand side is finite for 2 < 5 < 2*. Therefore, 


P-9 

V 




(3.3) 


f iV-2 

/ |a:| ^ 2 \vn\^dx ^ 0, as n —>■ 00 . 

JBi 


In order to obtain the same for the second integral in (j3.II) . we use the equivalence of the norm 
of Pir with the norm of 77,((M^), the radial subspace of 77^(M^), and the compact imbedding 


[ 

Jm.^ 


\x 




for any 2 < q <2*. Then, 



and finally 
(3.4) 


f |x| 2 dx < f |x| \Vn\‘^ dx ^ 0, 

B1 JbI 

as re —>• oo. ■ 

By unitary equivalence, the result holds in terms of re, hence the inclusion 

Hr ^ L9(M^), 2<q<2*, 

is compact, where Hr is the subspace of H consisting of radial functions. 

The following weighted interpolation inequality is immediate by setting a = /3 = j in 

the classical paper [7]. 

Lemma 3.2 Let 2 < q < and N >3. Then, there exists a constant C = C{N, q) > 0, such 
that 

N(q-2) 

(3.5)/" \x\~‘^^~\v\‘^dx < C ( f ) ( f \x\~^^~‘^^v\‘^dx 

Jm \Jm J 

for any v £ H. 

Remark 3.1 Notice that inequality liS. ,5|) holds without any hypothesis on radial symmetry. 


2q-N(q-2) 

) ‘ 


3.2 Global Solution 


We quote the basic points of the standard theory of semilinear Schrodinger equations (cf. [9l [10] ) 
ensuring the global well-posedness on the initial value problem 


iift + 1 ^ + IV’I'' V = 0 , 

- 0 ( 0 ) = - 00 . 


The operator defined by 
(3.7) 


D{C) := {^GHr-.Cy^e X}, 

Cip := |x|^-2v • (|x|-(^-2)v<^). 


defines a self-adjoint operator with T < 0. Therefore, by unitary equivalence, the operator defined 
by 


(3.8) 


DiC) := {'iP £ Hr : Cip £ X} , 

£0 := A0+(A^)"^, 


is also self-adjoint with £ < 0, where 0 = |x| Furthermore, iC defines a skew-adjoint 

operator and generates a group of isometries on Hr- Adapting step by step the arguments of m 
Theorem 3.3.9], we may easily conclude the following: 


Theorem 3.1 Let N > 3 and 2 < q < -^ + 2. Then, for all 0o G Hr, there exist T = r(0o) > 0 
and unique solution 


f;£C{[0,T),Hr)nC\[0,T),H-^), 
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of liS. 6\) . where ^ is the dual space of Hr- In addition, for all t £ [0, T), the following properties 
hold: 

(3.9) 


and 

(3.10) 

Set 


/ \'ij:{t)\‘^dx = / I'tpol'^dx (conservation of charge), 

E{'if{t)) = E^ipo) (conservation of energy). 


E{u) '■= - [ {ul'^dx. 


Interpreting inequality (13.5p in terms of u, it follows that there exists e G (0,1), such that 


(3.11) 


Fiu) < 


1 — £ 


WWh + C'(ll''^llL 2 (KiV)), for all u € Hr- 


Therefore, in analogy to [9l Theorem 3.4.1], we may set T{'tpQ) = oo in Theorem 13.11 for all 
V'o G Hr- 


3.3 Stability 

We introduce the functional J : Hr 


defined by 
1 , 


J{v) := E{\x\ 2 u) + 


|x|-(^-2)|u|2dx. 

Clearly, the minimization problem (|1.8p is equivalent to the following, 

(3.12) 
where 

~ r , , ri-2 

T ■.= {v £ Hr '■ \x\ 2 V £ r}. 

Both minimization problems are well-defined since, by inequality (|3.5p and the property 2 < q < 
^ + 2, there exist (5 > 0 and K < oo such that 


u G T, 

J{v) = min{J( 2 :) : z G T}, 


(3.13) 

Let 


J{v) > 5 ||u||^ — K, for all u G T. 


k-y := inf{J('(;) : v £ T}. 


In the following lemma we establish the precompactness of any minimizing sequence with respect 
to ky. 

Lemma 3.3 Let 7 > 0. Then, any sequence (u„) in Hr satisfying 
(3.14) J{Vn) —>■ ky, J'{Vn) “>■0, Vn £ T, 

contains a convergent subsequence in Hr- Furthermore, its limit solves the minimization problem 

(E2D. 
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Proof. First, observe that for n large enough, 


h.y -\- o(l) ^ JQ ^ ^ ^ P 


(3.15) 


g-2 

2q 


■^riWrir ’ 


that is, Il'ynll’H,, is bounded. Going if necessary to a subsequence, still denoted by (vn), there 
exists V € Tir such that ^ u in PLj.- By Lemma IHITI — >■ u in \x\~^^ 2 ~dx). 

Also, notice that 


= < J'{Vn) - J'{v),Vn - V > 


(3.16) 


+ / kl ‘^Vn-\v\'^ ‘^v){Vn-v)dx. 


Evidently, 

(3.17) < J'{vn) — J'{v),Vn — u > —)■ 0, as n —>■ oo. 

Moreover, it follows by Holder’s inequality that 


(3.18) 


[ |x| '?^2^(|?;„|9 ^v^-\v\ 1 ^v){Vn-v)dx <c( 

JRN 


x\ ^ 2 \y — y 




\ 1/9 


which tends to zero, as n —)> oo. Therefore, we found a u G PLr and a subsequence (un) such that 


(3.19) Vn^V, mPLr- 

It is obvious that u G F. Finally, by the weak lower semicontinuity of J and the definition of k.y, 
we obtain J{v) = k.y, and the proof is complete. ■ 

By a solution to (jl.2p we mean a couple {X.y,u.y) G M x 77^, where A..,, is the Lagrange multiplier 
associated to the critical point of E on F. If v.y is a global minimizer of J on F, then there 
exist {\.y,v.y) G M x Pir solving the elliptic equation 

(3.20) - V • + X.y\x\-^^-^^V.y - \x\-'^^\v.y\‘^-^V.y = 0. 

By unitary equivalence, S'.y / 0 and e*^'>'*|x| ^ v.y corresponds to a standing wave of (jl.ll) . 

Observe that if u.y G S^, then G S^, for all t > 0. Also, is a periodic function in 

time, hence we may say that S.y consists of a class of closed orbits. In this context, the following 
definition of (local) orbital stability makes sense. 


Definition 3.1 The set S.y is said to be orbitally stable, if for any e > 0, there exists a 6 > 0, 
such that for any global solution 'tp{t) of id. 6\) with 

(3.21) dist{'ipQ, S.y) < 6, 

it holds that 

(3.22) 

where 

(3.23) 


dist{if{t), Sy) < e, for all t > 0, 
dist{w,S^):= inj \\w — z\\jj^. 
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Proof of Theorem 11.11 Arguing by contradiction, assume that there exist sequences (V’On) C 

Hr, (tn) C M+, and eo > 0 with 

(3.24) llV’On ~ 

for some u (z S^, and such that the global solutions ifn with initial values ^fQn satisfy 


(3.25) 

Let us set 

Clearly, if 

then 


inf Wifnitn) - z\\jj > eo- 


zes. 


Un ■■= 'tpn{tn)- 


k-y := min{Li(u) : u € L, u is radial} 


— k-y 


We deduce from (I3.24p that 

(3.26) E{'ijjQn) ^ K and / \ijjon\'^dx ^-f. 

Jrn 

Applying the conservation laws (13.9p and (|3.1Up . we obtain 

(3.27) E{un) —>■ ^7 and / \un\‘^dx —>■ 7 . 

Jrn 

We choose a sequence /3„ ^ 1, such that 

(3.28) Ei^PnUrf) ^ ^7 and / dx — 7 ; 

JRN 

that is, fdnUn is a minimizing sequence of problem (II. 8 p . For example, we may insert 

\ 1/2 

7 \ 


(3.29) 


/3n = 


(_ 2 _) 

KIrJV \Un\^dx) 


Setting Vn = \x\ 2 Un, we obtain 

(3.30) J{PnVn)^k^ and / \x\~^^~'^'^\(3nVn\'^dx = j. 

JrN 

Ekeland’s variational principle [361 Theorem 2.4] yields to the existence of another sequence 
Cn G F satisfying 

(3.31) J{Cn) k^, J'iCn) 0, IlCn — PnVn\\-^ < — • 

nr 

By Lemma [ 3 ^ we deduce that there exist a subsequence in Hr, still denoted by (n, and a C G <S '7 

N-2 

such that Cn \x\ 2 (" in Hr- Therefore, 

inf \\Un-z\\H^ < \\Vn-\x\^C\\Hr 


(3.32) 


— I|W T PnVn Cn T Cn |^| ^ CIlHr 

1 N—2 

< |1 “/3n| llwll-^ H 1“ llCn — l^;] 2 ClI'Hr^O) 

nr yj 


which contadicts p3.25p and the proof is complete. ■ 
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3.4 Behavior at the Origin 

In this section, we describe the behavior of the standing wave solutions, obtained in Subsection 
E31 around the origin. 

Proof of Theorem \1.S\ Consider a global minimizer E of problem (11.81) for some 7 > 0, 
which can be chosen to be nonnegative on Then, there exists a Lagrange multiplier A-y, 

such that the pair (Ay,Vj) € Mx'Hj. solves the stationary problem (l3.2Up . where Vy = 

We introduce the following transformation 

~ 1 

(3.33) Wy(t) = Vy(r), t = (—logr) v -2 ^ r = \x\^ 

for 0 < r < 1. Note that 

lim w^(t) = lim v^(r). 

t^+00 r^l- 

Then, Wy satishes 

-Awy + X{N - 2)'^Vl{t)Wy + {N- 2)^V2{t)\Wy\'^-'^Wy = 0, 


where 

Vi(t)=exp(- 2 t-(^- 2 ))t- 2 (^-i) and = exp - 2 ) . 


If we set Tl(0) = 1 ^ 2 ( 0 ) = 0, then Vi, V 2 are continuous functions. Standard regularity results 
imply that Wy{0) is a well-defined real number, while Maximum Principle implies that this number 
cannot be zero. Thus ^-,,(0) is well-defined and positive, which means that Uy behaves at zero 
exactly as . ■ 

Consequently, the results of [33l Section 2.5] yield the presence of a correcting term in the energy 
norm. The exact value of the norm of Uy is given by the formula 

~ “ A.[Uy) -|- ||'U-),||^2(']gJV) , 


where 


A{Uy) = 


N{N - 2) 


UN Vy{0). 


On the contrary, we shall examine in next section a case where the new energy term acts in an 
additive way to the total energy. 


3.5 The Case of Hardy Energy at Infinity 

We consider the following Cauchy problem 

. 334 ^ f i\x\-^wt + Aw+{^f = 

\ w{0) = Wo, 


where r(; = r(;(|x|), N > 3 and 2 < q < jj + 2. The approach is based on the unitary equivalence 
with problem (13.6p . To this end, we introduce the Kelvin transformation in the form 

'ipiy) = \x\^-'^w{x), * = 


(3.35) 
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denoted hy ip = K.{w). Clearly, for smooth functions, it holds that 


Ay'ijj{y) = \x\^^^Axw{x) and 


'tpiy) 


= ixr+2 


w{x) 


|y|" \x\ 


and the equations are equivalent for y ^ 0. The differences appear in the energy near the 
singularity versus the energy at infinity. 

As it was mentioned formerly, the energy space corresponding to (|3.6p is Hr with norm given 
by (ini). We shall use this formulation for the definition of the energy space W corresponding 
to (I3.34p . Indeed, following exactly the proposal stated in [33], the space W is defined as the 
isometric space of Hr under transformation K, and the correct energy norm is given by the 
following formula 

llrcll^ = lim(/i/j(tc) + Ai/£-(tc)) + / \x\~^\w\'^dx, 

efO jRiV 

where ^ 

h/eiw) := [ \Vw\^dx- f ^Adx, 

V 2 J |x|2 


1/e 


and 


^i/e{w) := ^ [ \w\^dS. 

2 J\x\=l 


More precisely, W is defined as the completion of the C'^(M^\{0})-functions under the norm 
(3.36) llrcll^ =/rjvPu;) + / \x\~'^\'w\'^dx. 


Remark 3.2 The question of well-posedness, as well as, the existence and orbital stability of 
standing waves for \3.34^ in the spaee W is understood through the unitary equivalenee with Hr- 


Remark 3.3 Note that in the case of a function w, which behaves at infinity like |x| ^ 

we have again the appearance of a correction energy term with measure 

(3.37) AM = 

where v{x/\x\^) = |x|“2“ w{x). However, there is a notable difference with problem 113.6]) . since 
in this case the singularity effect acts in an additive way to the usual Hardy functional. 


4 The General Case 


In this section, we are concerned with the orbital stability of the Schrodinger operator, equipped 
with the critical inverse square potential, without assuming any symmetry hypothesis. We con¬ 
sider the problem 


iipt + A'ip+ 1^ + g{x)\'4)\^ V = 0, 

V'(O) = V’o, 


where g G is a nonnegative function. 

The approach is an adaptation of the method used in the radial case. Therefore, we make an 
outline of the basic steps. 
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w > -iV + at 0, 

00 < —N + ^^^2 


Lemma 4.1 Let l<q<2*. If g satisfies 
(4.2) g r‘^, with 

then, the inclusion H ^ Lg{M.^, |a;| ’i~^dx) is compact, where 

|x| 2 “'^dx) := G L^(M^) : J |f I'^dx < oo| . 

Proof. Let Vn be a bounded sequence of C'^(M^\{0})- functions in H; without loss of generality, 
we assume that Vn converges weakly to 0. Then, we split the norm of Vn as 

WVnWu = \\Vn\\n{Bi) + I I'^^n | |'H(B=) • 

In addition, we study separately the radial parts = v^(r), r = |x|, and the nonradial parts 
u”'’ of Vn- The proof consists of four steps. 

Step 1: The radial part uj) in Bi. Arguing as in lemma ItTI we obtain that u)) G II^(Bi C M^), 
whenever G TL^Bi). Therefore, 


f |x| \v‘^\^ dx ^ 0, for any pG(l,+oo), 

Jbi 

as n ^ oo. Then, for p large enough, 

f I I rwj f I n r, 

/ |a:| 2 = \x\ ' ' g\\x\ p dx 

Jb^ Jb, 

\x\-^^-^Uv:\pdx] . 


< 


jj;! ^ ' ^(P-P^ gp-q dX 


'Bl 


'Si 


Let g ~ at 0. Then, the hrst integral of the right hand side is hnite if 


(4.3) 

Finally, we get that 


a; > (iV - 2)| - iV, 


f Al-2 

/ |x| 2 'I g\vl^y dx ^ 0, as n —>■ oo, 

Jbi 

for g( ~ r'^, at 0, and oo satisfying (|4.3p . 

Step 2: The nonradial part vffi in Bi. For the nonradial parts vlfi' of Vn, we observe that 
(4.4) \x\-^^-^'>/^ vff G 

This follows as in [35] or as in [211 Pg-196]; setting u'lfi' = , we have 


I nr||2 _ T 
\Ur, W'lJ — -tro 


'Jn WV. — + [ 

However, using decomposition into spherical harmonics. 




/ |V^ 


u'Zl'^dx, 


O'! a. 
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and ()4.4p holds. Moreover, 

so is also bounded in Thus, G H^{Bi), with H^{Bi) to be compactly 

embedded into LP{Bi), 1 < p < 2*. 

For e > 0 small enough, we set 
(4.5) 


Aq,e = 11 - 


2 * - e 

It is clear that is increasing as a function of e, so 


Aq,e > (1 - — 


-1 


for any e > 0 small enough. Then, 


L 


x\-^'ig\v';,^\Ux = [ g\\x\ 

Bi Jbi 


N—2 

■ — v^^\Ux 


< 


iBi 


g dx 


1 

^q,e 


'Bi 


Q 

I nr\2*-e ' 

\x\ ^ '^n ' 


dx 


Let 5 ~ r‘^ at 0. Then, the first integral of the right hand side is finite if ()4.3p holds. Finally, we 
get that 


/ 

JBi 


\x\ 2 g\v'^^\^ dx ^ f), as n —)■ oo. 


for g( ~ r‘^, at 0, and uj satisfying (14.3p . 

Step 3: The radial part in B^. In the case of the exterior domain B^, both in the radial and 
the nonradial case, we use the following Kelvin transform. We set 


(4.6) 


X 


v{x) = w{y), y = TT 2 


The determinant of the Jacobian of the Kelvin transformation in dimension d > 2 is equal to 
— \x\^^ and 

2 


^xw ; ,^2 

xr 


= kl Wyw{y)Y 


Then, 

(4.7) 

iBi JBi 

We restrict now in the radial case. For the radial functions w"^, we have 


^\\niBf)= [ \y\ f \y\ ^ ^\wn\'^dy. 

J J 


'Bi JBi 


Wn\ dy > ||w^nll//l(BiCR2)- 


\y\ p converges to 0 in L^{Bi), p > 1. 


Hence, as in step 1, 
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Moreover, 


f |x| 2 1 gdx = 

Jb9 


iBi 




iBi 


^tzln-OMj-E^n iV-2 

|yM ^ 9\\y\ ^ Kl'^dy 


< 




iBi 


iBi 


ixr(^-2) \v^jpdx 


9^ 

q 

p 


Let g ~ 1x1“^ at oo. Then, the first integral of the right hand side is finite if 


or, for p large enough, 
(4.8) 

Finally, we conclude 


-w > -q{N - 2)^^ + fV + -fV, 

2p p 




[ \x\-^'ig\v^JUx^0, as n ^ oo, 
Jbi 


for g ^ r‘^, at 0, and ui satisfying (|4.8p . 

Step 4- The nonradial part in B^. Using transformation (|4.6I) . we have that 
(4.9) \\vn\\H= [ iVwnpdy + / \y\~^~^\wn\'^ dy. 

Moreover, as in Step 2, we may obtain 

\\Vn\\n> f \^i\yr‘^<lfdy+ [ \y\-*U-‘^wZ\^dy. 

JRJV JrJV 

Restricting ourselves in Bi, we have 

|y| ~ w"^ is bounded in H^{Bi), 


so 


\y\ 2 w'^ converges to 0 in LP{Bi), l<p<2*. 


Then, 


[ \x\-^'i g\vll^\Ux = [ |y|V 9 - 2 iV^|y^nr|g 

Jb 9 Jbi 


dy 


[ \y\ 

JBi 




< 


iBi 


'Bi 


i\y\^N-2)g-2N 




Q 

2*-e 
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where is defined in (14.5p . Let g ~ \x\^ at oo. Then, the first integral of the right hand side 
is finite if 

N N -2 

-u>2N-{N -2)q-^> N - —q, 


A 


q,e 


or Lj satisfies 


Finally, we get that 

N-2 

\^\ ^ ^ dx ^ 0^ as n ^ oo, 


/ 

is? 


for g ^ r‘^, at oo, and uj satisfying ()4.8p 

Remark 4.1 (i) Since 
(4.10) 


N-2 

— N -\ -- —q < 0, 


for every 1 < q < 2*, we conclude that g maybe constant at the origin. On the other hand, it 
should decay to zero at infinity. As we saw in Lemma HOI this is not the case if we restrict on 
radial functions. 

(a) Consider the functions u = v. Then, 


I 


lx| 2 'ig\v\'^dx= / g\u\'^dx. 

Jrn 

In the special case q = 2, g should be a function such that 

g ^ r‘^, with 

as it is natural from the Hardy’s inequality. 


oj > —2, at 0, 
uj < —2, at oo. 


(Hi) Any function g with g € D L^*-i (M^) satisfies 

By Lemma 13.21 and the modified behavior of g at zero, the following is immediate: 
Corollary 4.1 Let 2 < q <2*, N >3. If g satisfies 

/ N ,, f w > 0, at 0, 

(4-11) , with < . 

[ + at oo, 

then, there exists a constant C = C{N,q) > 0, such that 


[ 

Jr^ 


g\x\ ^ 2 \v\^dx < 


\-{N-2) 


\Vv^dx 


N(q-2) 

4 


(4.12) 

for all V UTL. 


<^-^^\v?dx 


2q-N(q-2) 

4 


The energy functional naturally associated to (|4.ip is defined by 

■= - A£(^^)) - ^ [ gH'^dx. 

Z £ 4.0 q 
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Theorem 4.1 Let N > 3 and 2 < q < -^ + 2. If g satisies then, for all ipo G H, there 

exists a unique solution 


Ip G C{[0,oo), H) n C^{[0,oo), H ^), 

of (EZP, where H ^ is the dual space of H. In addition, for all t >0, the following conservation 
laws hold: 

(4.13) f \ip{t)fdx = f \^po\^dx, 
and 

(4.14) Egim) = Egi^Po). 


Consider the minimization problem 


(4.15) 


f u e r, 

\ Eg{u) = mm.{Eg{z) : z G T} , 


which is well-defined due to Corollary 14.11 By using similar arguments to those used in Lemma 
13.31 we may obtain the existence of a standing wave solution G LI of (14.ip . for a given 

7 > 0. If 

*^9.7 ■= {u G H : u is a solution of (14.150 } 

denotes the set of global minimizers, then a nonradial version of Theorem 11.11 is stated in the 
following: 


Theorem 4.2 Assume that N > 3, 2 < q < -^ + 2 and g satisfies condition 1^4.11 ). 
set Sg^r^ is orbitally stable for problem dp. 


Then, the 


4.1 Improved Hardy-Sobolev Inequality 

As already said in the Introduction, there is no possibility for a Hardy-Sobolev inequality to hold. 
This is clear since functions behaving at the origin like , do not belong to L^*. As in 

the bounded domain case, we have to consider Improved Hardy-Sobolev (IHS) inequalities; i.e. 
we have to find a weight function h{x) such that 

/ r \{N-2)/N 

(4.16) \\(P\\h>{ I h{x)\(pf] , for any </> G C'^(M^). 

Unfortunately, in the case of the whole space we cannot apply the method described in [34p Lemma 
9.1]. This method provides a way of constructing optimal IHS inequalities (i.e. it provides h{x)) 
and it was used in [32l [34l |39] . However, we are in position to present an IHS inequality, although 
it is not optimal. Before this, we remark the following: 

1. The nonradial part of cp belongs to (see (j4.4jl L so in this case (I4.16[) holds with 

h{x) = 1. 

2. The radial part (pr of (p satisfies ()4.16p with h[r) = jxp. This may be obtained by following 
the arguments of Lemma l3.H in this case \x\^^cpr G LL^(ffi^). 
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3. For any R> 0 and (p G inequality (14.161) holds with 

2(N-1) 

h{x) = log ^ , D> R. 

This optimal (IHS) inequality was proved in [2T]. In the case of radial functions, D is 
allowed to be equal to R (see [39]). Note that, since D or R appear in the inequalities, 
these inequalities depend on the domain. 


Directly from 1. and 2. we have the following. 


Lemma 4.2 Assume that h is defined as 


h{x) 


|a;p, lx| < 1, 

1, |x| > 1. 


Then, inequality is true. 


However, from the point of view of 3. this inequality is not optimal. 
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